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Using a semi-classical approach, we describe an on-chip cooling protocol for a micro-mechanical 
resonator by employing a superconducting flux qubit. A Lorentz force, generated by the passive 
back-action of the resonator's displacement, can cool down the thermal motion of the mechanical 
resonator by applying an appropriate microwave drive to the qubit. We show that this on-chip cool- 
ing protocol, with well-controlled cooling power and a tunable response time of passive back-action, 
can be highly efficient. With feasible experimental parameters, the effective mode temperature of a 
resonator could be cooled down by several orders of magnitude. 

PACS numbers: 85.85.+j, 45.80.+r, 85.25.Dq 

I. INTRODUCTION 

The rapid development of nano-technology has enabled fabrication of micro or nano-mechanical resonator^ with 
high frequency, low dissipation and small mass. In the way to observe quantized mechanical motion-, thermal 
fluctuation has become a major obstacles. The limited cooling efficiency and poor heat conduction at the milli-Kelvin 
temperatures of cryogenic refrigerators has stimulated a number of studies on the active cooling of micromechanical 
resonator (MR) in both classical regime and quantum regimes^^^^ ' 10 ' 11 ' 12 ' 13 ' 14 ' 15 ' 16 ' 17 ' 18 ' 19 ' 20 ' 21 . Of these proposals, 
the optomechanical cooling is the most highly developed one experimentally. For example, it has been successfully used 
to demonstrate cooling^ 9 - from room temperature down to 6.82 mK and from 2.2 K down to 2.9 mK. Optomechanical 
cooling has also been studied theoretically in the quantum limi t 22 i 23 i 24 . 

In a thermal environment, a MR randomly vibrates around its equilibrium position due to thermal fluctuation. 
The effective temperature of a MR mode can be defined by the mean kinetic energy of this mode. Cooling of MR 
is equivalent to suppressing the mean amplitude of the Brownian motion, which is the major barrier to precise 
displacement measurement. The optomechanical cooling of an MR is achieved by the passive back-action or active 
feedback of such optical forces as bolometric force 7 -^ and radiation pressure — i^iL^uLIii from a laser driven optical 
cavity. 

Instead of employing an optical component, in this paper, we present a cooling protocol for an MR in an on-chip 
superconducting circuit containing three or four Josephson junctions (JJ). The three-junction superconducting loop 
forms a flux qubi t 26 i 27 ' 28 . The typical energy scale of the Larmor frequency of the flux qubit is normally several GHz 
which is much larger than the oscillation frequency of the MR studied in this paper. The persistent current in this 
loop exerts a Lorentz force on an MR under an in-plane magnetic fiel d 29 ' 30 . This circulating superconducting current 
is modulated by the thermal motion of the MR through nonlinear Josephson inductance in a delayed way. Thus the 
Lorentz force exerted on the MR depends on the motion of the MR itself. This force acts as a passive back-action on 
the MR as the radiation pressure in the optomechanical cooling strategy. The thermal motion of the MR is damped 
by this back-action force with appropriate parameters. To use a microwave bias to drive the superconducting circuit 
helps to take away the thermal energy of the resonator. In other words, in the linear response regime, a microwave 
drive together with a flux qubit in a dissipative environment can be treated as an effective bath whose temperature 
is lower than that of the environmental temperature To. Thereby the mode temperature of the resonator, which is 
proportional to the mean kinetic energy of the fundamental oscillation mode, is decrease d 14 ' 31 ' 32 . The strong and 
tunable coupling between superconducting current and the MR displacement is favorable as regards obtaining highly- 
efficient and well-controlled cooling. Based on feasible experimental parameters, we also provide a detailed analysis 
of the cooling efficiency of this scheme. The estimation shows that it could constitute a promising alternative to 
optomechanical cooling. 

This paper is organized as follows: In Sec. II, we describe the setup we use to implement our cooling scheme. 
Sec. Ill presents the intuitive picture of this back-action cooling by Lorentz force, together with a general formalism 
for dealing with self back-action cooling. In Sec. IV, the cooling efficiency of this physical system is determined 
based on a detailed calculation of the "spring constant" and the response time of the Lorentz force. The lengthy 
calculation part concerning to the master equation and its steady-state solution are given in the Appendix. In Sec. V, 
cooling efficiency is estimated by using feasible experimental parameters . Finally in Sec. VI, we discuss the possible 
advantages, fluctuation and measurement protocol of this cooling protocol. 
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FIG. 1: (Color on line) Schematic diagram of our setup. A doubly-clamped mechanical beam is incorporated in a supercon- 
ducting loop with a superconducting Josephson junction (indicated with crosses) flux qubit. The initial bias in the loop is 
controlled by magnetic flux $ z o in z direction. A coupling magnetic field Bo is imposed on the beam in the y direction and the 
supercurrent under this magnetic field imposes a Lorentz force on the beam. A microwave line introduces a microwave bias to 
the qubit loop. The magnitude of the Lorentz force depends on the motion of the beam through the change of the supercurrent 
in the loop in a delayed way. This delayed back-action damps the motion of the beam in the z direction and thereby cools 
down the thermal motion of the beam. 

II. THE SETUP 

A schematic diagram of the setup for this cooling protocol is shown in FigfT] In the x-y plane, a doubly-clamped 
micro-mechanical beam with an effective length Lq is incorporated in a superconducting loop with three small- 
capacitance Josephson junctions. This mechanical beam can be created from a surface-micromachined silicon or 
carbon nanotube^ coated with superconducting material, or a self-supporting metallic airbridge. The fundamental 
vibration mode of the beam can be well approximated by using a harmonic resonator with oscillation frequency tub- 
With a proper bias magnetic flux, two classical stable states of the 3-JJ loop carry persistent currents in opposite 
directions. There is a finite tunneling rate A between the two classical persistent-current states (throughout this article, 
we let h — 1). By choosing parameters carefully, the subspace spanned by the two states is well separated from other- 
energy levels and the superconducting loop with Josephson junctions forms a two-level system (superconducting flux 
qubit) 26 ' 27 and coherent dynamics can be observe d 34 ! 35 . Recently, the microwave-induced cooling of a flux qubit has 
been demonstrated experimentally 36 assisted by the third energy eigen-state of the three-junction superconducting 
loop. The qubit ground state \g) and excited state |e) are coherent superpositions of two persistent current states 
denoted by |0) (clockwise current state) and |1) (anti-clockwise current state). The energy spacing between the two 
eigenstates is Q — VA 2 + e 2 where e = 2I P (Q e xt — $o/2) is the energy spacing of the two classical current states (|0) 
and |1)), with <fr ex t the external magnetic flux through the loop, I p the largest persistent current in the loop and $ 
the flux quantum. A microwave line is placed close to the circuit and generates microwave drive with frequency u>d 
on the flux qubit. Under the coupling magnetic field B a along y direction, the persistent superconducting current 
generates a Lorentz force Fl on the MR along z direction. This force couples the flux qubit with the oscillation 
motion of the MR. In the quantum regime of the MR, this configuration provides a solid-state analog of cavity QED 
system in strong coupling limit 29 . Similar setting with a coupled MR and dc-SQUID was proposed recently to study 
the displacement detection and decoherence of mechanical motio n 37 ! 38 ' 39 

In the present paper, we concentrate on a different regime where the oscillation frequency of the MR is much smaller 
than the Larmor frequency of the flux qubit so that the MR can be treated as a classical harmonic oscillator. Since 
the qubit and the MR energy scales differs greatly, the dynamics of the composite system can be handled following 
the line of the Born-Oppenheimer approximation: The master equation of the qubit is established by assuming a 
certain MR displacement z; the steady-state solution for the qubit dynamics is inserted into the classical Langevin 
equation of the MR to derive the noise spectrum of the mechanical displacement. 

III. THE COOLING MECHANISM 

In this section, we first present an intuitive understanding of the cooling protocol based on the Lorentz force back- 
action. Using the classical Langevin equation of the MR, we then proceed to analyze how this back-action leads to 
the suppression of the random motion of the MR. The explicit form of the Lorentz force back-action in the Langevin 
equation will be derived in the next section. 

The composite flux qubit and MR system is kept at an environmental temperature Tq in a dilution refrigerator. The 
initial constant bias flux of the loop is Q z o in the z direction. When the MR is displaced to z(t) by thermal noise, the 
total magnetic flux in the superconducting loop is changed to $ 2 o + BoLgzit). In principle, for an MR with an aspect 
ratio close to 1, thermal fluctuation also induces oscillation in the y direction. However, since the magnetic field in 
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the z direction is much smaller than the coupling magnetic field Bq in the y direction, in this work we only study the 
motion in the z direction. An increase or decrease in magnetic flux leads to a change in the persistent current / in 
the loop after the system reaches a metastable state after a response time delay r rosp . Owing to the nonlinearity of 
Josephson junction inductance, the supercurrent in a metastable state depends on the total bias magnetic flux and 
hence depends on the displacement of the MR. Since the Lorentz force is proportional to the supercurrent, the Lorentz 
force exerted on the MR Fl(£ + T reS p) depends on the displacement of the MR z(t) before the time delay. This means 
that there exists a passive back-action mechanism for the MR: the motion of the MR leads to a delayed force on 
the MR itself. We found that, with a proper microwave drive (red-detuned with the qubit), this passive back-action 
damps the thermal motion of the MR. This cooling protocol is similar to that of the self cooling experiments based 
on optomechanical coupling 7 ' 25 . 

The above intuitive picture can be clearly understood by studying the MR dynamics. To accomplish this, we start 
by establishing an equation of motion for the MR. 

The coupling term in the Hamiltonian of qubit-resonator composite system is 2 ^ 

#i„t = B L Iz. (1) 
where / is the current operator. Then the Lorentz force on the resonator is 

F L = = - W (2) 

oz 

Suppose the MR is displaced at z, after a time delay T rcsp , the resulting Lorentz force in the z metastable state 
(z) = Tr (^pq S ^ (z) Fj^j depends on z, where is the reduced density matrix of the qubit in a metastable state. 

For a small displacement z, the Lorentz force can be expanded to the linear order of z as (z) = Fq + k^z with 
the effective "spring constant" of the Lorentz force. Since this force is a delayed response to the motion of the 
MR, its effect on the MR can be described by a delayed response function h (t — t') = 1 — e~ 7 ( t-t ' with 7 = l/r reS p. 
Under the Lorentz force, the equation of motion for the MR mode with mass to, rigidity k = mui 2 and an inherent 
damping rate V is as follows^ 

d 2 z dz , f* dF { r s) [z (*')] n , / 

where F t h is the Brownian fluctuation force which is related to the environmental temperature (F t h (t) F t h (£')) = 
2/cbTotoF(5 (t — t'). The stable solution of the motion equation j3|) is related to the effective mode temperature T e g 
of the MR by the equipartition theorem k e g (z 2 ) = /csT c ff where ks is the Boltzman constant. Then the cooling 
efficiency 77 = Tq /T e g can be written as 



rc c 2 ff r+™ dLo 



K 2 ff -^ 2 ) 2 + r 2 ff u 



(4) 



where 

F cff = r + ri, 



7 2 



0J 2 + 7 2 k 



(5) 



and 



r rn _Ubl_k L . , 

ui 1 + "f z k 



with Qm = Wft/r is the quality factor of the MR. For \ki,\ <C k, the integral can be explicitly carried out as£ 

k u^+Y 

From Eq.([7]), it can be seen that the sign of fcz, determines whether the resonator is cooled or heated, i.e. when k^ 
is positive, the Lorentz force damps the motion of the MR and the final temperature is lower than the original one, 
and vice versa. For a positive fc/,, the cooling efficiency increases linearly with k^. However, it should be noticed that 
both the linear response regime and the definition of a single-mode resonator are valid only for khjk <C 1. 
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The cooling efficiency is shown by a contour plot in Fig. [2] as a function of two ratios Wb/7 and k^/k. Wb/7 
characterizes the ratio between the time scale of the response time and the oscillation period while k^/k is the scaled 
cooling strength of the passive Lorentz force. The magnitude of the efficiency is indicated by the gray level: Higher 
efficiency is represented with a lighter color. It is seen that the cooling efficiency increases with k^/k. The strong 
coupling in the solid state cavity QED makes it promising to achieve larger k^. For a given k^/k, the largest cooling 
efficiency can be achieved by optimizing cjb/7. The optimal points for each k^jk value are indicated by the red dots 
in Fig. O The optimal uJb/j slightly increases with k^. For k^/k <C 1, which is in the case of interest, the red dots 
indicate that the optimal cooling is realized for cob = 7. This means the largest cooling efficiency is achieved when 
the back-action response time matches the oscillation period of the resonator. As we show later, the response time 
is of the order of the relaxation time of the flux qubit. The qubit relaxation rate ranges from sub MHz to several 
tens of MHz, depending on the operating point e u 40 i 41 , this implies that the flux qubit loop is suitable for cooling a 
resonator over broad frequency range. Moreover, since the operating point £0 can be controlled by magnetic flux <fr z o, 
the back-action response time is tunable in situ. 



IV. CALCULATION OF THE BACK-ACTION OF THE LORENTZ FORCE 



As we discussed in the previous section, without considering additional fluctuation, the damping effect of the Lorentz 
force leads to cooling or heating depending on the sign of the effective spring constant k^\ A positive k^ leads to 
cooling while a negative k^ results in heating. The cooling efficiency is proportional to the magnitude of k^. Before 
proceeding with our discussion of the cooling efficiency in realistic experiment, it is necessary to obtain an explicit 
expression of kh from the dynamics of the flux qubit system. 

The Lorentz force in a metastable state for a given displacement z is 

F { L S) (z) = -B I P L (a z ) s (8) 

where (cr z ) s is the expectation value of a z in a steady-state. And a z and a x are Pauli matrixes defined by the persistent 
current states as: 

a z = |0)(0|-|l)(l| 

^ = 10x11 + 11x01 (9) 

(<r z ) can be computed from the dynamics of the driven flux qubit in a dissipative environment. To do this, we start 
with the effective Hamiltonian of the qubit at a given displacement z of the MR 

e (z) A 

H q = -^TT-Vz + ~^<*x + A(Jz C0SUJ dt- (10) 




FIG. 2: (Color on line) The dependence of the cooling efficiency on the ratios of fci/fc and Wb/7. The magnitude of the cooling 
efficiency is re-scaled and indicated by the gray level (higher efficiencies are shown lighter) . The other parameters are the same 
as those used for estimating the cooling efficiency in Sec. V. The red dots indicates the point with the largest efficiency at each 
given fct/fc value. 
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Here A characterizes the amplitude of the microwave drive, 

e(z) = e + gz (11) 

with g = 2BqI p Lq, Eq = 2/ p <I> z o denotes the initial bias away from the degeneracy point. If the qubit is biased at the 
degeneracy point, £o = 0. 

By defining a new set of Pauli operators 

a' z = a z cos 9 (z) + a x sin 9 (z) , 

a' x — —a z sin 9 (z) + cr x cos 9 (z) , (12) 
we diagonalize the first two terms of eq. (TlT)]) as, 

H q = ^y^-<t' z + A {a' z cos9{z) - a' x sin6> (z)) cosuj d t, (13) 

where 

cos0( 2 ) = g|, sin^).^, (14) 

and 



fi(z) = ^A 2 +£ 2 (z)taSl +gz, (15) 



with g = geo/rio, Qq = \/ A 2 + e§ is the energy spacing between the qubit eigenstates when the displacement z of the 
MR is zero. Here we have used the fact that the displacement z is very small and cosujdt is a fast-oscillating term. 

If the drive is near-resonant with the qubit frequency £7q, by performing unitary transformation Ur = eyLp(ia' z LOdt /2) , 
the Hamiltonian (|13[) in the rotating frame is transformed to 

Hr = i (|^(*)) Kit) + U R (t) HUUt) = ^4 + ^< (16) 

with A' (z) = — Asin9 (z), 5 (z) = Suj + gz and i5oj = i^o — is the detuning between the qubit free energy and 
the drive. The term Aa' z cos 9 (z) cos ojdt and Asm9(z)a + e luJdt + h.c in eq. (Tj"3"|) is neglected because exp(±iu>dt) and 
exp(±2iw^) are fast-oscillating terms since the drive frequency is close to the qubit energy spacing. The diagonalized 
Hr 



Hn = (17) 



is expressed by another set of Pauli operators 



with 



and 



a z = <j' z cos (3 (z) + a' x sin (3 (z) , 

a x = —a' z sin (3 (z) + a' x cos (3 (z) , (18) 



sm(3(z) = ^f\, cos/3(z) = ^, (19) 



uo (z) = v/<5 2 (z) + A' 2 (z) . (20) 

With the above transformation relations, we can derive the master equation for this dissipative flux qubit in a 
rotating frame. Solving the master equation to obtain its steady-state solution and using the above relations to 
transform it back to an experimental frame, we can obtain (a z ) s in an experimental frame (see Appendix for detail) 



(21) 
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From the last equation of ([38)1 , the response time needed for the supercurrent to reach this value is 

1 _ 1 2 

7 ~~ T 1 + cos 2 /?, 



T rcsp — _ — ^ , _j a (22) 



with cos/3o = cos/3(z = 0). The response time is of the order of the qubit relaxation time. Shifting the qubit bias will 
change the qubit relaxation time and hence change the response time. On the other hand, as found with SET-assisted 
cooling^, the response time also depends on the drive detuning and power. This is in contrast to optomechanical 
cooling where the response time is always determined by the cavity ring-down time. 

Since z is small, the force can be expanded to the linear order of z: (z) = F + ki,z, where 

F = -B I P L Q 2 ^^ cos O , (23) 
1 + cos^ po 

and the effective spring constant of the Lorentz force is 

k L = B I p L (^^) • (24) 
The derivation of the above equation can be explicitly carried out as 



z=0 



_ (2B I P L ) 2 ( el 2uj SujA 2 uj A 2 Suj 2 \ 



u \ (wg + Suj 2 ) 2 ^0 ^0 w o + Slj 

where loq is the value of eq. (|20|) at z = 



^,= J ( W--A 2 ^. (26) 



When the flux qubit is biased near the degeneracy point, e < A, in the bracket of eq. 1(25]). the second term is much 
smaller than the first one owing to the small prefactor uio/flo (in this case, Slo is about several gigahertz while the 
detuning and drive amplitude are about several megahertz) and we neglect the second term 



k L « (2B I p L f ^ 2 - ( 27 ) 

"o (wq + Sur) 

Eq. (|25p represents the explicit form of expressed by the physical quantities of this system. As we have already 
noticed, the cooling efficiency is closely related to this effective spring constant k^. Using this expression, we are able 
to study the cooling efficiency of this physical system. Inserting eqs. (|25|) and (|22|) into eq. 0, we obtain the cooling 
efficiency of this system 

_ l | (2B I P L ) 2 4Q M r (1 + cos 2 go) / £ 2 2lo 5ujA 2 | lu A 2 5lo 2 \ 
71 mujbivo 4w 2 + T 2 (1 + cos 2 /3 ) 2 \^1{u} 2 +Suj 2 ) 2 f!^ 2 + fo 2 j' 

r) — 1, i.e., the second term of the right part of the above equation, is the change in the cooling efficiency caused by 
the Lorentz force back-action. If 77 — 1 is positive, the back-action leads to cooling, and vice versa. The dependence 
of 77 — 1 with respect to the drive detuning 5u> and amplitude A near degeneracy point is shown in Fig. [3] It can be 
seen that for a positive Suj (red detuning), the cooling regime is reached, and vice versa (however, note that 77 < 
corresponds to an unstable case rather than heating). There is no cooling or heating for zero drive ^4 = 0. This is 
consistent with optomechanical cooling. However, because of the difference between SU(2) and Heisenberg algebra, 
for Suj > 0, the cooling efficiency does not scale with A monotonically but reaches its maximum at A — \/2Suj as shown 
in Fig. |4j Physically, this is because too strong a drive destroys the back-action mechanism owing to its dominant 
rapid dynamics and even drive the qubit out of the two- level subspace. It is also shown in Fig. [3] that, as A increases, 
the variation in the cooling efficiency with respect to Suj around the peak tends to become much flatter. This feature 
makes the cooling efficiency robust as regards operating point fluctuation by increasing the drive amplitude. However, 
the largest cooling efficiency also decreases with the increases in A. Therefore, in real experiment, a trade-off is needed 
to obtain both stable control and high cooling efficiency. 

It is worth to mentioning that, if the qubit is biased at the degeneracy point, the second term in eq. (|25p cannot 
be omitted. Thus we obtain a small but finite cooling effect even at the degeneracy point. Based on the first-order 



2SujA 2 
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perturbation theory, side-band cooling is not possible at the qubit degeneracy point. However, driven Rabi oscillation 
at degeneracy point has been demonstrated experimentally 4 - . This is explained by fluctuation of the operating point 
or higher-order perturbation theory 43 . It is interesting to observe that this effect is preserved by our semi-classical 
treatment of back-action cooling. However, based on the treatment of quantized MR, other resonant conditions may 
be needed if we are to study this phenomenon^ 3 -. 



V. COOLING EFFICIENCY IN PRACTICAL SYSTEMS 



Having obtained the effective spring constant of the Lorentz force and established the relationship between the 
cooling effect and the back-action mechanism, we can now estimate the cooling efficiency of this protocol for a real 
system. To achieve higher cooling efficiency, it is necessary to increase the mechanical quality factor Qm and the 
coupling strength k^. A higher quality factor implies that the MR is less affected by the thermal environment. Greater 
coupling strength means stronger back-action damping. Moreover, as shown in Fig. [21 another crucial factor with 
respect to improving the cooling efficiency is the ratio between the MR frequency and the flux qubit relaxation rate, 
i.e., the time scales of the oscillation period and the feedback response from the qubit. Efficient cooling is achieved 
when the two time scales match. In self-cooling experiments with optomechanical coupling, this is one of the main 
experimental challeng esii In our proposal, the response time depends on the relaxation rate of the flux qubit. This 
rate can be controlled from 0.1 MHz to 10 MHz by controlling the initial bias magnetic flux <I>zo in the z direction. By 
controlling the high frequency ( A/H) noise, e.g., by deliberately attaching an external impedance to the circuit, we 
can change the qubit spontaneous emission time. Hence the relaxation rate of the flux qubit can be further monitored 
to match the frequency of the mechanical mode. This feature enables us to optimize the cooling of MR with different 
mode frequencies and greatly improves the cooling power. 

We made our estimation based on a flux qubit loop with an aluminum Josephson junction. In this case, the 
relaxation time around the degeneracy point of the qubit is of the order of one hundred nanoseconds to several 
microsecond a 40 ' 41 . Hence it is suitable for a cooling MR mode with MHz frequency. For example, at an environmental 
temperature of 300 mK, a flux qubit with I p = 600 nA, A = 5 GHz, Eq — 1 GHz and the relaxation rate close to the 
degeneracy point is 7 = 5 MHz. The qubit is driven by a microwave with u>d — 5.089 GHz and A = 14.1 MHz. When 
the driven qubit is used to cool a doubly-clamped Si beams^ with cj& = 5 MHz, k = 0.1 N/m, Lq = 5 /mi and quality 
factor Qm — 10 4 , assuming Bq — 5 mT, we obtain fc^ = 2.72 x 10 -3 N/m. The cooling efficiency r\ is about 1.27 x 10 2 , 
which means the effective mode temperature of the MR can be cooled to about 2.36 mK. Since k^ is proportional to 
the square of the coupling magnetic field, increasing Bq can greatly enhance the cooling power. In principle, the upper 
limit of Bq is the critical magnetic field B c of the superconducting material in order to preserve superconductivity. 
For aluminum, B c « 9.9 mT. But the in-plane magnetic field can be much stronger than this critical value (e.g. even 
larger than 100 mT— ) without destroying the superconductivity. In addition, suitable arrangement of the magnetic 
field can result in the strong magnetic effect acting only in the MR region while being largely canceled out at the 
junction location. This means this cooling protocol could be potentially more powerful than the above estimation. 




FIG. 3: (Color online) The dependence of the change in cooling efficiency rj — 1 on the drive frequency 8lo for different drive 
amplitude. Suj and A are in milli-Kelvin units. All the other parameters are the same as those used to estimate the cooling 
efficiency in Sec. V. 
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FIG. 4: (Color online) The dependence of the change in the cooling efficiency 77 — 1 on the drive amplitude A with drive detuning 
5uj = 1 mK and 5u> — 5 mK respectively. 5u> and A are in milli-Kelvin units. All the other parameters are the same as those 
used to estimate of cooling efficiency in Sec. V. 

However, one major problem with too strong an in-plane magnetic field in the experiment is the possible tiny 
vibration of the sample with respect to the coupling magnetic field. With a strong magnetic field, a significant 
shift may be induced in the operating point by the tiny vibration. One way to solve this problem is to design the 
gradiometer-type flux qubit or to use an on-chip magnetic field generating coil 2 ^. 



The basic idea behind our cooling protocol is similar to the self-cooling of a micro-mirror by optomechanical 
coupling 44 : The flux qubit plays the role of the optical cavity and the microwave drive on the qubit acts as the 
laser drive on the cavity. The Lorentz force produces a passive back-action on the resonator in the same way as 
photothermal force or radiation pressure. 

However, the mathematical treatment of a two-level artificial atom is rather different from that of a bosonic field. 
Moreover, in practice the two systems have distinct physical nature: (1) An on-chip solid-state system without any 
optical component might have certain advantages as regards its application; (2) The coupling between a flux qubit 
and a resonator can be controlled by controlling the applied magnetic field Bq which is independent of the qubit free 
Hamiltonian and microwave drive. This is different from the coupling system of a Cooper pair box and NAMR 47 
where the bias voltage V g modifies the coupling coefficient as well as the free Hamiltonian of the charge qubit. This 
feature gives the system more flexibility in terms of increasing the cooling power and switching the cooling process on 
and off ; (3) In order to near-resonant to laser, the frequency of the FP cavity is more than 1 THz. To achieve efficient 
optomechanical cooling for a MHz oscillator, the optical quality factor of the FP cavity should be about 10 s . This 
calls for a mirror with extremely high finesse. But in our case, it is easier to match the two energy scales since the 
relaxation rate at the degenerate point is typically of the order of microsecond. (4) The relaxation rate of a flux qubit, 
which determines the back-action response time, can be modified in situ by the bias magnetic flux in the z direction. 
In principle with a GHz oscillator, we can work in both the 7 « u>b and 7<wt regimes for this Lorentz-force cooling. 
The latter regime is desirable for cooling towards the quantum limit 22 . This implies that our present proposal might 
be able to cool the MR into its quantum ground state. However, this question requires further discussion involving a 
consideration of quantum fluctuation. 

In the above discussion, we only considered the damping of the resonator caused by the qubit via the Lorentz 
force. However, according to the fluctuation-dissipation theorem, fluctuation is always associated with dissipation. 
Additional fluctuation is also induced by the back-action of the qubit system. For a more comprehensive description, 
we need to take this fluctuation into account. The additional fluctuation can be represented by an effective temperature 
T\. Therefore, the resonator is effectively in contact with two different bath a 14 ' 31 > 32 , one is the real environment with 
temperature To, and damping rate T and the other is the effective bath with effective temperature T\, and damping 
rate T%. Physically, the effective bath is formed by the lossy qubit under microwave drive. The temperature of the 
MR reaches a balance between the two baths as 



VI. DISCUSSIONS 



= TTp + riTi 
r + Tr 



(29) 
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The near-resonant microwave drive largely suppresses the fluctuation (see Appendix) 



Therefore we obtain T% ~ u>b <C Tq ~ f2, which means the fluctuation caused by the qubit can be effectively disregarded, 
as far as the spontaneous emission of the qubit is concerned. A more realistic estimation of this fluctuation requires 
detailed information about the quantum noise of the qubit. 

The cooling of the MR could be measured by using the conventional motion transduction methoo—,. In this system, 
it can also be revealed by the reduction in the integration of the power spectrum for the MR (z 2 (u))). Since 

2V2ASm B I p L ; 

01 = -7= — -=— =■ f — oz, (31) 

A 2 + 2Suj 2 A x J 

the motion power spectrum is proportional to the current power spectrum. The MR power spectrum can be recorded 
by detecting the current in the loop. The current in the flux qubit loop can be recorded with a dc SQUID switching 
measurement or with a more sophisticated phase sensitive dispersive readout such as Josephson bifurcation measure- 
ment^. 
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APPENDIX: THE DISSIPATIVE DYNAMICS OF THE QUBIT 

The dissipative dynamics of the qubit is governed by the following master equation 



p=-~[H q ,p]+£p (32) 



where p is the reduced density matrix of the qubit, H q is the qubit Hamiltonian in eq. (|10p and C is the Liouvillian 
characterizing the influence of the environment. 

As regards the qubit, we only consider the spontaneous emission with rate Tq and neglect the excitation and 
dephasing terms because the qubit energy spacing is much larger than the environment temperature and the qubit is 
biased close to the degeneracy point 

To 

Cp = — {2<j-p<7 + — p<T + <J- — o + <j-p) . (33) 

By successively applying three unitary transformations to eq. (|32|) 50 or making use of the Fermi golden rule^, we can 
readily obtain the master equation in the interaction picture of the rotating frame 

pf = £ R pf, (34) 

where 

£ R P? = y (2*-P? <?+ - pf <r+<r- - v+v-pf) 
{2d+pf a- - pf a-a+ - or_*+p?) 
+^-(a z pfa z -pf), (35) 
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with 

r x = ^(l + cos/?) 2 , 

T T = (1 - cos/3) 2 , 

T v = ysin 2 /3 (36) 

and To is the spontaneous emission rate. This master equation shows that the effective bath formed by the driven 
qubit in dissipative environment is 

exp (^UEi = p + cos M 2 07) 

Therefore, in the interaction picture the equations of motion for the average value of the qubit operators are 

| (>- W > - - gt±a±°a (>-(«)>, 

I ft (<)) = - (r t + r, ) (s, (<» - (r, - r, ) . (38) 

Setting d(<Ji)/dt = 0, we can obtain the steady-state solution of the qubit in a dissipative system 



(t+) s = (*->. = 0, 

r T -r ; 2 cos /3 (z) 



r T + r x i + cos 2 ^(z)' 



(39) 



Making use of eq. (|39|) and the transformation relations eqs. (fT2|) . (fT6|) and l|18[). it turns out that the quantity of 
interest in the experimental frame is related to those in the rotating frame as 

(<r z ) s = {cos9(z) cos /3 (z) — sin# (z) sin/3 (z) cosuidt) (<5z) s ■ (40) 

Neglecting the high-frequency oscillating term, we obtain the effective value in eq. (l2~Tj) 

2cos 2 /3(z) 



o- z „ = - 



1 + cos 2 /3 (z) 



cosO(z) (41) 
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